Effects of Strain on Orbital Ordering and Magnetism at the Perovskite Oxide 

Interfaces: LaMnOg/SrMnOg 



00 

o 
o 

(N 



o\ 



C/3 



C^ 



O 

o 



> 

o 

oo 
O 



X 



B. R. K. Nanda and Sashi Satpathy 
Department of Physics & Astronomy, University of Missouri, Columbia, MO 65211 

(Dated: July 9, 2008) 

We study how strain affects orbital ordering and magnetism at the interface between SrMnOa 
and LaMnOa from density-functional calculations and interpret the basic results in terms of a three- 
site Mn-O-Mn model. Magnetic interaction between the Mn atoms is governed by a competition 
between the antiferromagnetic superexchange of the Mn t2g core spins and the ferromagnetic double 
exchange of the itinerant eg electrons. While the core electrons are relatively unaffected by the 
strain, the orbital character of the itinerant electron is strongly affected, which in turn causes a 
large change in the strength of the ferromagnetic double exchange. The epitaxial strain produces 
the tetragonal distortion of the MnOe octahedron, splitting the Mn-e^ states into x'^-y^ and 3z^-l 
states, with the former being lower in energy, if the strain is tensile in the plane, and opposite 
if the strain is compressive. For the case of the tensile strain, the resulting higher occupancy of 
the x^-y'^ orbital enhances the in-plane ferromagnetic double exchange owing to the larger electron 
hopping in the plane, causing at the same time a reduction of the out-of-plane double exchange. 
This reduction is large enough to be overcome by antiferromagnetic superexchange, which wins to 
produce a net antiferromagnetic interaction between the out-of-plane Mn atoms. For the case of the 
in-plane compressive strain, the reverse happens, viz., that the higher occupancy of the 3z^-l orbital 
results in the out-of-plane ferromagnetic interaction, while the in-plane magnetic interaction remains 
antiferromagnetic. Concrete density-functional results are presented for the (LaMn03)i/(SrMn03)i 
and (LaMn03)i/(SrMn03)3 superlattices for various strain conditions. 

PACS numbers: 75.70.Cn, 71.20.-b, 73.20.-r, 71.70.-d 



I. INTRODUCTION 

Recent advances in successfully designing atomically 
sharp interfaces between dissimilar transition metal ox- 
ides have revealed the formation of new electronic and 
magnetic phases at the vicinity of the interface, which 
are qualitatively different from the parent compounds. 
The interfacial phases show diverse magnetic properties 
due to the coupling between charge, orbital and spin de- 
grees of freedom. For example the magnetic ordering at 
the interface between the two antiferromagnetic insula- 
tors SrMnOs (G-type) and LaMnOs (A-type), schemat- 
ically shown in Fig. [TJ could be ferromagnetic along all 
directions, ferromagnetic in the xy-plane and antiferro- 
magnetic normal to the plane, or antiferromagnetic in the 
plane and ferromagnetic normal to the plane depending 
on the composition of the parent compounds and epitax- 
ial strain on the interfacei'^'^'^'^i^. 

The epitaxial strain, arising due to lattice mismatch 
between the constituent compounds of the superlattice 
and the substrate, induces anisotropic hopping between 
orbitals to cause orbital ordering at the interface. By 
varying the strain condition the orbital ordering changes 
which in turn changes the magnetic ordering at the inter- 
face. In this paper we examine the magnetic properties 
at the interface of SrMnOa (SMO) and LaMnOg (LMO) 
for different epitaxial strain conditions through first prin- 
ciples electronic structure calculations. 

Experimental studies show that if the substrate in- 
duces tensile strain at the interface of the LMO/SMO 
superlattice, where the in-plane lattice parameter 'a' is 



greater than the out-of-plane lattice parameter 'c', as 
in the case of (LMO)3/(SMO)2 superlattice grown on 
SrTiOa (STO) substrate, the magnetic ordering of the in- 
terfacial Mn atoms is A-type with in-plane (Mn02 plane) 
ferromagnetic (FM) ordering and out-of-plane (between 
Mn02 planes) antiferromagnetic (AFM) ordering^. 

Quite interestingly, when the (LMO)3/(SMO)2 super- 
lattice is grown on Lao.3Sro.7Alo.65Tao.35O3 (LSAT) sub- 
strate, which induces no strain (a ~ c), the interface 
shows a three dimensional FM ordering (F-type)i. If 
the interface experiences a compressive strain (a < c) , as 
in the case of LMO/SMO superlattice grown on LaAlOs 
(LAO) substrate, the magnetic ordering is C-type with 
in-plane AFM ordering and out-of-plane FM ordering^. 

Substrates are instrumental in inducing epitaxial strain 
and thereby enforce tetragonal distortion to the super- 
lattice. As a consequence, in case of LMO/SMO su- 
perlattice, the substrate distorts the MnOg octahedron 
and splits the degenerate Mn-e^ states into x^-y^ and 
3z^-l states. Varied tetragonal distortion changes the 
onsite energy and hence the occupancy of these two non- 
degenerate eg states (Fig. [5]). Since the electronic config- 
uration of Mn atoms away from the interface is the same 
as in the bulk compounds, (Mn*+, tLe") for SMO and 



n3+ 



tLcg) for LMO, strain is not expected to affect 



the magnetic configuration of the inner Mn02 layers to 
a large extent. However, at the interface, where we see 
the valence state of the Mn atoms lies between 3-1- and 
4-1- because of charge reconstruction^!^, the varied occu- 
pancy of the non-degenerate e^ orbitals imposes different 
orbital ordering for different strain conditions (Fig. [2|) 
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FIG. 1: Different magnetic configurations considered in tliis 
paper for tiie (LMO)i/(SMO)i superlattice. 'A' stands for 
tlie structure witli ferromagnetic ordering in the Mn02 plane 
and antiferromagnetic ordering between the planes, while 'F' 
stands for ferromagnetic ordering in all directions and 'C 
stands for antiferromagnetic ordering in the Mn02 plane and 
ferromagnetic ordering between the planes. The schematic or- 
bital ordering shown in the figure was found from our density 
functional results presented below and was also inferred from 
the experiments^. The symbols J and J' denote, respectively, 
the out-of-plane and in-plane exchange interactions between 
the Mn atoms. The strain condition under which each struc- 
ture is stabilized has been indicated in the figure. The oxygen 
atoms which occur at the mid-point between two neighboring 
Mn atoms have not been shown. 
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FIG. 2: Energy splitting of the Mn(eg) orbitals at the 
LMO/SMO interface for compressive and tensile strain condi- 
tions. The parameter A is the difference between the energies 
of the Mn-d^2_y2 and the Mn-d322_i orbitals. Compressive 
strain makes the 3z^-l orbital lower in energy while tensile 
strain makes it higher. 



and influences the interface magnetism considerably. 

In this paper, we have studied in detail the interfa- 
cial magnetic properties of LMO/SMO superlattices for 
different strain conditions by performing electronic struc- 
ture calculations based on the density functional theory. 
To illustrate the strain effect on magnetism, we have pro- 
posed a simple three-site model to calculate the intcr- 
facial Mn-O-Mn magnetic exchange both in the Mn02 



plane and between the planes for different strain con- 
ditions. From the model we see that the onsite energy 
difference between x^-y^ and 3z^-l orbitals (Fig. [J) is 
instrumental in switching the ferromagnetic and antifer- 
romagnetic interactions. When the 3z^-I orbital is suffi- 
ciently lower in energy than the x^-y^ orbital (compres- 
sive strain), the Mn-O-Mn exchange is antiferromagnetic 
in the plane and ferromagnetic between the planes, and 
opposite when x^-y^ orbital is sufhciently lower in en- 
ergy (tensile strain) . If the energy levels of both the e^ 
orbitals are close enough (lattice matched interface) , then 
the Mn-O-Mn exchange is ferromagnetic in all directions. 
The rest of the paper is organized as follows. In sec- 
tion II we describe the structural and computational de- 
tails. A detailed analysis of the electronic structure of the 
(LMO)i/(SMO)i superlattice at different strain condi- 
tions, obtained from the density- functional calculations, 
is carried out in section III. In section IV, we illustrate the 
effect of epitaxial strain on the magnetic ordering, with 
the aid of a proposed three site (Mn-O-Mn) model. Elec- 
tronic and magnetic properties of the (LM0)i/(SM0)3 
superlattice at different strain conditions are discussed in 
section V. Finally in section VI we present the summary. 



II. STRUCTURAL AND COMPUTATIONAL 
DETAILS 

We have taken the equivalent cubic perovskite struc- 
ture of LMO and SMO in order to study the electronic 
and magnetic properties at the interface of these two 
compounds with the aid of first principles electronic 
structure calculations. The effect of epitaxial strain, 
which arises due to lattice mismatch between the sub- 
strate and the LMO/SMO superlattice, is taken into ac- 
count by applying tetragonal distortion to the superlat- 
tice. 

The tetragonal distortion is quantified by the 'c/a' ra- 
tio which differs from one. Here 'a' is the in-plane (xy- 
plane) lattice parameter which coincides with the lattice 
parameter of the substrate and 'c' is the average out-of- 
plane lattice parameter (along z-axis). The 'c/a' ratio 
is determined from the linear relation: c - ao = -4z^(a - 
ao), where 'ap' is the in-plane lattice parameter of the 
superlattice when there is no strain (c/a = 1) and coefff- 
cient V is the Poisson's ratio which is approximately 0.3 
for perovskite manganites^'^. Experimentally it is found 
that for LMO/SMO superlattices, ao matches with the 
weighed average of the lattice constants of bulk LMO 
(3.936 A) and bulk SMO (3.806 A)l. For example, for 
(LMO)i/(SMO)i superlattice, the value of ao is i(3.936 
+ 3.806) A. 

In this paper, we have considered two superlattices, 
viz., (LMO)i/(SMO)i and (LM0)i/(SM0)3 to study the 
electronic and magnetic properties at different strain con- 
ditions. As is well-known, the strength of the Jahn- Teller 
distortion is less in the mixed compounds (La, Sr)Mn03 
as compared to that of LaMnOa, we have considered a 



small Jahn- Teller distortion (Q2 « 0.05 A) in the basal 
plane for the interfacial Mn02 layers. However, test cal- 
culations showed that a small variation of Q2 does not 
change the electronic and magnetic properties of the su- 
perlattice qualitatively. 

All electronic structure calculations reported in this 
work have been performed using the self-consistent 
tight-binding linearized muffin-tin orbital (TB-LMTO) 
method with the atomic sphere approximations (ASA)^. 
Self-consistent calculations are done within the frame- 
work of generalized gradient approximation including 
Coulomb correction (GGA+U). All results arc obtained 
with U = 5 eV and J = 1 eV unless otherwise stated. 



III. ELECTRONIC STRUCTURE OF THE 

(LaMnOs ) 1 /{SrMnOs ) 1 SUPERL ATTICE 

In this section, we describe the effect of strain on 
the electronic structure at the interface from ab ini- 
tio density-functional (DFT) calculations. We focus on 
the (LMO)i/(SMO)i superlattice and our results sug- 
gest that many of the interfacial electronic and magnetic 
properties shown by this superlattice should also be valid 
for the more general (LMO)„/(SMO)m superlattices. 

We briefly summarize the electronic structure and 
magnetism for the bulk SMO and LMO compounds. In 
bulk SMO, the Mn atoms are in A+ charge state to have 
three d-electrons which are occupied in the triply degen- 
erate t2g states. The doubly degenerate Cg states, which 
are higher in energy with respect to t2g states because 
of a MnOg octahedral crystal field split, remain unoccu- 
pied. The tjo spin majority states mediate an antiferro- 
magnetic superexchange to stabilize the G-type antifer- 
romagnetic ordering in the bulk SMO compound, where 
spin of each Mn atom is opposite to that of the nearest 
neighbor Mn atomsiSiiii. 

In bulk LMO the Mn atoms are in 3+ charge state 
to have four d-electrons. Three electrons are occupied 
in the localized t2g states and the remaining one elec- 
tron is occupied in the e^ state. The Jahn- Teller (JT) 
distortion to the MnOg octahedron further splits the Bg 
states into two non-degenerate states: e^ which is lower 
in energy and ei which is higher in energjsi^. The one 
eg electron is occupied in the e^ state whose lobes are 
pointed towards the longest Mn-0 bond. The JT distor- 
tion stabilizes the A-type antiferromagnetic structure in 
the LMO compound^*. 

At the LMO/SMO interface the Mn atoms do not sat- 
isfy the 4+ charge state or the 3+ charge state to support 
the bulk magnetism of SMO or LMO. The mixed valence 
nature of the Mn atoms as well as the effect of epitax- 
ial strain create diverse magnetic phases at the interface, 
which will be analyzed in this section. 

Epitaxial strain, arising due to the substrate on which 
the interface is grown, induces tetragonal distortion to 
the cubic interface which is quantified by the 'c/a' ra- 
tio that differs from one. Experimental studies show 
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FIG. 3: Total energies for the magnetic configurations, A 
and C, relative to the energy for the configuration F, as a 
function of the tetragonal distortion c/a (top). The magnetic 
configurations A, F, and C are shown in Fig. [l] Bottom figure 
shows the occupancies of the x^-y^ and the 3z^-l orbitals per 
Mn atom as a function of the tetragonal distortion. 



different magnetic behavior at the interface for different 
'c/a' ratios^. To obtain the dependence of the magnetic 
ground state on the strain condition, we have performed 
total energy calculations in the range 0.95 < c/a < 1.05 
for three possible magnetic configurations (A, F, and C) 
(Fig. [1]) . Magnetic configuration A represents the FM or- 
dering in the Mn02 plane and AFM ordering between the 
planes. Magnetic configuration C represents the AFM or- 
dering in the Mn02 plane and FM ordering between the 
planes and F represents the FM ordering in all directions. 
The energetics are shown in Fig. [3] (top). 

From the figure we see that for a strong compressive 
strain {e.g. c/a — 0.95), 'A' is the most stable magnetic 
configuration. For the lattice matched structure ( c/a 
= 1, no strain), the interface stabilizes with magnetic 
configuration F and in case of a strong compressive strain 
(e.g. c/a = 1.05) it stabilizes with magnetic configuration 
C. The results are in accordance with the experimental 
observations which show that when the substrates are 
STO (c/a = 0.98), LSAT (c/a = 1.01), and LAO (c/a 
= 1.05), the respective magnetic configurations at the 
LMO/SMO interface are A, F, and C.-i 



We see that as strain changes, the occupancy of the e^ 
orbitals, which controls the magnetic interaction at the 
interface, also changes. This is shown in Fig. [3] (bottom) . 
For the tensile strain condition (c/a < 1) the occupancy 
of the x^-y^ orbital is greater than the occupancy of the 
3z^-l orbital and opposite if the strain is compressive (c/a 
> 1). For the lattice matched structure (c/a — 1) both 
the Bg orbitals are more or less equally occupied. Fig. [3] 
also shows that for any value of 'c/a', the non-degenerate 
eg states combinedly occupy 0.5 electrons which along 
with three t2g core electrons make the average valence of 
the interface Mn atoms to be +3.5 as expected. 

Magnetic interaction between the Mn atoms is deter- 
mined by the competition between ferromagnetic double 
exchange- '^'^-''^- via the itinerant Mn-Cg electrons and an- 
tiferromagnetic superexchange between the localized Mn- 
t2g core spins. When x^-y^ is more occupied and 3z^-l 
orbital is less occupied (or unoccupied) , the strong double 
exchange in the Mn02 plane strengthens the ferromag- 
netic ordering while superexchange stabilizes the antifer- 
romagnetic ordering between the planes. The magnetic 
ordering is opposite to the above when the occupancies 
of the two Cg orbitals are reversed. If both the e^ orbitals 
are more or less equally occupied, the double exchange 
stabilizes the ferromagnetic ordering both in the plane 
and between the planes. 

As described in the following subsections, a detailed 
analysis of the density-functional electronic structure of 
the LMO/SMO interface under different strain conditions 
gives us a better understanding on the strain induced 
orbital ordering and its effect on magnetic properties at 
the interface. 



A. c/a = 0.95, Tensile strain 

Tensile strain reduces the out-of-plane lattice parame- 
ter c and enhances the in-plane lattice parameter a. In 
other words it decreases the Mn-0 bond length between 
the Mn02 planes and increases it in the plane. In such 
a scenario, the total energy calculation (Fig. [3]) suggests 
a stable A-type magnetic configuration (Fig. [1]) when 
tetragonal distortion (c/a) is close to 0.95. In Fig. 21 
we have shown the total and partial densities of states 
(DOS) for the (LMO)i/(SMO)i superlattice (c/a = 0.95) 
in the A-type structure obtained from the GGA-fU cal- 
culations. 

The characteristic features of the electronic structure 
under tensile strain as seen from Fig. [5] are as follows. 
The localized Mn-t2g states lie far below the Fermi level 
(Eip) because of the octahedral crystal field and strong 
Coulomb repulsion. The 0-p states occur in the energy 
range of -6 to -1 eV. The x^-y^ and 3z^-l orbitals are 
predominant at E^^. Since the intraplane (on the xy- 
plane) Mn-0 bond is longer than the interplane (along 
the z-axis) one (Fig. [2]), this lowers the energy of the x^- 
y^ orbital making it more occupied and raises the energy 
of the 3z^-l orbital, which becomes less occupied. 
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FIG. 4: Total and partial DOS for the (LMO)i/(SMO)i su- 
perlattice (c/a — 0.95) in the A-type magnetic configuration 
The symbols t and i represent the local spin of the atoms 



The Mn-Bgl state at the Fermi level (Ef) splits into x^ "^ 



and 3z -1 states. The orbital character of the e^ states at Ef 
is shown in Fig. [5] 



The origin behind the stability of A-type magnetic con- 
figuration for the tensile interface is explained below. In 
the bulk LMO, Mn (3-I-) atom has the electronic con- 
figuration tlg^l and in the bulk SMO, Mn (4-I-) atom 
has the electronic configuration tjgCg. Since at the in- 
terface, the Mn02 layers are surrounded by (SrO)" layer 
and (LaO)^"*" layer, the interface Mn atoms arc left with 
the average valence state of +3.5. In such a scenario, the 
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FIG. 5: (color online) . Orbital character of the electron bands 
near Ef for the three magnetic structures A, F, and C. The 
bands are plotted along the high symmetry points F (0, 0, 0), 



X(^, 



0), Z (0, 0, -^),andR(^, 



, -j^). The unit 
cell for the magnetic structures is doubled along the xy-plane 
with the formula unit 2 x (LMO)i/(SMO)i. For the AFM 
configurations (A and C), we have two Mnf and two Mn| 
atoms. Only spin-up bands are shown for the FM configura- 
tion (F). 
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FIG. 6: (color online). Valence electron charge- density con- 
tours plotted on the xz- plane in the energy range Ef - 0.15 
eV to Ef to indicate the orbital ordering for the A-type mag- 
netic configuration. Contour values are: Pn ~ po x 10" e/A^ , 
where po = 3.7 x 10~^, S = 0.4 and n labels the contours. The 
charge contours on the yz-plane (not shown) are identical to 
that of xz- plane. The orbital ordering is mainly x'^-y'^. 



t2g orbitals will occupy 3 electrons in the spin majority 
states and the e^ orbitals will occupy the remaining 0.5 
electrons. 

Without any occupancy of the e^ states, the only con- 
tribution to the energy comes from the superexchange 
interaction between the localized t2g states to stabilize 
the G-type AFM phase as in the case of SMO. However, 
the itinerant eg states, if partially occupied, can mediate 
the Anderson-Hasegawa double exchange^^ to stabilize 
the FM phase. 

The strength of the FM ordering in the plane or out 
of the plane depends on the occupancy of the individual 
x^-y^ and 3z^-l orbitals. From our calculations (Fig. [3]) 
we find that for tensile strain condition (c/a — 0.95), 



the occupancy of x^-y^ orbital is close to 0.45, while for 
3z^-l orbital it is less than 0.1. This is also reflected 
in the Mn-e^ band dispersion for the A-type magnetic 
configuration shown in Fig. \E\ Since, the unit cell for the 
magnetic structure is doubled along the xy-plane (i.e. 2x 
(LMO)i/(SMO)i), for the AFM magnetic configuration 
we have two Mnt and two Mn| atoms. Hence, for the 
local spin majority channel, we have two x^-y^ orbitals 
and two 3z^-l orbitals. From the figure we see that the 
3z--l orbitals are mostly in the conduction band and only 
one x^-y^ orbital of two crosses the Fermi level and lies 
mostly in the valence band. This implies that almost 
one electron per two Mn atoms in the x^-y^ states is 
occupied which is consistent with the orbital occupancy 
calculation. 

In such a case the x^-y^ orbitals will mediate the dou- 
ble exchange mechanism in the Mn02 plane to stabilize a 
FM ordering in the plane. The gain in kinetic energy due 
to the planar orbital order, induced by the anisotropic 
hopping, is more than the loss of super exchange energy. 
Since the 3z^-l orbitals are only marginally occupied, su- 
perexchange between the localized t2g electrons stabilizes 
the AFM ordering between the Mn02 planes. The net 
result is an A-type AFM ordering at the interface. 

The valence electron charge density contours for states 
in the vicinity of the Fermi level (Ef), shown in Fig. [SI 
provides a visualization of the above analysis. The charge 
contours show that the orbital ordering is predominantly 
Mn-x^-y^, 0-pa; and Py, while the occupancies of the 3z^- 
1 and P2 orbitals are small. As a result we see a strong 
coupling between the Mn-e^ and 0-p orbitals in the plane 
while it is rather weak between the planes. Therefore, 
the in-plane magnetic exchange interaction J' is ferro- 
magnetic while the out-of-plane J is antiferromagnetic 
(Fig. [T]). Our results are consistent with the experimen- 
tal results that the magnetic ordering at the interface for 
(LMO)3/(SMO)2 supcrlattice grown on STO substrate 
(c/a = 0.98) is A-typei. 



B. c/a=1.0, Lattice matched structure 

Lattice matched interfaces are without any tetragonal 
distortion and hence the in-plane and out-of-plane Mn- 
O bond lengths are identical. Total energy calculation 
(Fig. [3]) in this case favors a three dimensional FM or- 
dering (F-type). To gain insight into the origin behind 
the FM ground state, we analyze the electronic structure 
for the lattice matched interface. In Fig. [7l we have 
shown the total spin-up and spin-down DOS for the F- 
type magnetic configuration. 

General features of the electronic structure of the lat- 
tice matched interface are similar to that of the tensile 
interface. However, now on either side of the the Fermi 
level, both x^-y^ and 3z^-l orbitals are predominant in 
the spin- up channel and they have nearly equal onsitc en- 
ergies. It is due to the fact that the Mn-0 bond lengths 
are same both in-plane and out-of-plane, making the e^ 
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FIG. 7: Total spin-up and spin-down DOS for 

(LMO)i/(SMO)i superlattice in the F-type magnetic 
configuration. Both x'^-y^ and Sz'^-f orbitals are more or less 
equally occupied. The orbital character of the eg states at Ef 
in the spin-up channel is shown in Fig. [5] The superlattice 
shows the half-metallic behavior. 
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FIG. 8: (color online). Valence electron charge- density con- 
tours plotted on the xz-plane in the energy range Ef - 0.15 
eV to Ef to indicate the orbital ordering for the F-type mag- 
netic configuration. Contour values are: Pn = po x 10" e/A'^, 
where po = 3.7 x 10~ , S — 0.4 and n labels the contours. The 
charge contours on the yz-plane (not shown) are identical to 
that of xz-plane. The orbital state of each Mn is a mixture of 



-y2 and Sz^- 



1. 



states nearly degenerate in energy. This is also substan- 
tiated from the dispersion of the spin-up Mn-eg bands for 
the F-type structure shown in Fig. [5l Since the formula 
unit is doubled along the xy-plane, there are four Mn 
atoms and all are in the same spin orientation. Hence, 
in the spin-up channel, we have eight eg bands of which 
almost six lie in the conduction band. Of the remaining 
two bands, which are part of the valence bands, one is 
predominantly of 3z^-l character, while the other is pre- 
dominantly x^-y^. Hence the occupancy of each of these 
orbitals is close to a quarter electron per Mn atom which 
is also seen from the orbital occupancy results of Fig. [3l 



The valence charge-density contours of Fig. [8] indicates 
the orbital occupancy of the two Mn-e^ orbitals as well 
as their hybridization with the 0-p orbitals. 

The partially occupied x^-y^ and 3z^-l orbitals medi- 
ate a ferromagnetic double exchange, strong enough to 
overcome the antiferromagnetic superexchange both in 
the plane and out of the plane to stabilize a three dimen- 
sional FM ordering. We have shown earlier^^ that in the 
case of CaMnOs/CaRuOs interface, a leaking of 0.2 elec- 
trons from the metallic CaRuOs side to the Mn-Cg states 
near the interface, which were otherwise unoccupied, is 
sufhcient to stabilize the FM ordering of the Mn spins. 
In the present case, both the e^ orbitals being occupied 
substantially (more than 0.2 electrons in each orbital), a 
strong ferromagnetic double exchange coupling along all 
directions is expected. This is consistent with the experi- 
mental observation of ferromagnetism in the LMO/SMO 
interface structures grown on the LSAT substrate (c/a = 
1.01)1. 

The other prominent feature in the electronic structure 
of the lattice matched interface is the opening of a gap 
at the Fermi level in the spin-down channel which makes 
the system half- metallic (Fig. [7]). 



C. c/a = 1.05, Compressive strain 

When the strain is compressive, the Mn-0 bond length 
reduces in the Mn02 plane while it increases between the 
planes. As a result, the 3z^-l orbital is lower in energy 
and is more occupied, while the x^-y^ orbital is higher in 
energy and is less occupied which is seen from the densi- 
ties of states (Fig. ^ as well as from the band structure 
(Fig. [5l right panel). 

As in the case of A-type magnetic configuration in the 
tensile strain condition discussed earlier, here also we 
have two Mnf and two MnJ, atoms. So for the local spin 
majority channel, we have two x^-y^ orbitals and two 
3z^-l orbitals. From Fig. [5]we see that the x^-y^ orbitals 
lie in the conduction band and only one of the two 3z^-l 
orbitals lies in the valence band. This shows that the 
occupancy of the 3z^-l orbital per Mn atom is close to 
0.5 and x^-y^ orbitals are basically unoccupied. This is 
seen from the orbital occupancy (Fig. [3]) as well as from 
the charge-density contour plot of Fig. (TUl where we see 
that the orbital ordering is predominantly Mn- 3z^-l and 
0-p,. 

The partially occupied 3z^-l orbital mediates a strong 
double exchange mechanism to make the out-of-plane 
magnetic ordering ferromagnetic. The in-plane remains 
antiferomagnetic due to the superexchange between the 
localized t2g electrons. Thus the net magnetic configura- 
tion for the compressive interface is C-type in agreement 
with our total energy calculations shown in Fig. [3] Ex- 
perimental studies on (LMO)3/(SMO)2 interface, grown 
on LAO substrate (c/a — 1.05), do indeed show a C-type 
antiferromagnetic configuration^ consistent with our the- 
oretical results. 
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FIG. 9: Total and partial DOS for the (LMO)i/(SMO)i su- 
perlattice with compressive strain (c/a = 1.05) in the C-type 
magnetic configuration. The symbols | and J, represent the 
spin majority and minority states with respect to a Mn atom. 
The spin-majority eg band splits into x^-y^ (unoccupied) and 
3z^-l (partially occupied) bands as seen more clearly from the 
band structure plot (Fig. [S|. 
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FIG. 10: (color online). Valence electron charge- density con- 
tours plotted on the xz-plane in the energy range Ef - 0.15 
eV to Ef indicating the orbital ordering for the C-type mag- 
netic configuration. Contour values are: pn = po x 10" e/A'^, 
where po = 3.7 x 10~^, 5 = 0.4 and n labels the contours. The 
charge contours on the yz-plane (not shown) are identical to 
those plotted on the xz-plane. The orbital ordering is mainly 
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Unlike the lattice matched interfaces which are metal- 
lic, the compressive interface is insulating. It is known 
that the strong correlation effect in manganites plays an 
important role to drive the insulating behavior. To elu- 
cidate this effect, in Fig. [TT]we have plotted the Mn-d 
DOS for different values of U. For [/ > 5 eV, a gap opens 
at the Fermi level to make the interface insulating. 



FIG. 11: Atom projected Mn-d DOS for a single Mn atom in 
the spin-majority channel for the (LMO)i/(SMO)i superlat- 
tice with the compressive strain condition (c/a = 1.05) and 
the C-type magnetic configuration. The results are obtained 
from GGA-I-U with different values of U. For U > 5 eV the 
DOS has a gap at the Fermi level. The z^ — 1 state is split due 
to the interaction between the Mn atoms on adjacent planes, 
while the x'^ — y^ state is not, because of the AFM alignment 
of the neighboring Mn atoms in the plane. 



IV. STRAIN EFFECT ON MAGNETISM: 
THREE SITE MODEL 



So far, from the electronic structure calculations for 
the (LMO)i/(SMO)i superlattice under various strain 
conditions, we found that strain changes the relative oc- 
cupancy of the two Mn-Cg orbitals, which in turn affects 
the magnetic ordering in the structure. We found that 
the ordering is antiferromagnetic in the plane and fer- 
romagnetic out of the plane, if the strain is compres- 
sive and opposite, if the strain is tensile, while for the 
lattice-matched interface, the magnetic ordering is ferro- 
magnetic in all directions. In this section we develop a 
simple three site model consisting of the Mn-O-Mn atoms 
to further understand the effect of the strain-controlled 
orbital occupancy on the magnetic interactions at the 
interface. 

In bulk perovskite manganites, the tj. core spins in- 
teract via the antiferromagnetic superexchange. In ad- 
dition to this, the eg electrons mediate the ferromag- 
netic Anderson-Hasegawa double exchange between the 
core spins, which competes with the antiferromagnetic 
superexchange. The strength of the double exchange de- 
pends on which of the Cg orbitals is occupied because of 
the anisotropic hopping and these are modeled in this 
Section. 

In the LMO/SMO superlattice, as seen from the 
density-functional results, the Mn atoms at the inter- 
face are left with one extra electron (0.5 electrons per 
Mn) which occupies the itinerant e^ states. The epitax- 
ial strain splits the degenerate Mn-e^ states into x^-y'^ 
and 3z^-l states. From the model below, we will see that 
depending on the strain condition the relative occupan- 
cies of the interface x^-y^ and 3z^-l orbitals could switch 
a ferromagnetic interaction into an antiferromagnetic one 
and vice versa. 

In our model we consider the t2g electrons as classical 
core spins which are fixed at the Mn sites with negli- 
gible intersite hopping as compared to the itinerant Cg 
electrons, as has been used by many authors in the liter- 
ature. Hence the Mn-O-Mn double exchange is due to the 
hopping between the itinerant x^-y^ and 3z^-l electrons 
and 0-p electrons (Fig. [T^ . The model Hamiltonian 
thus reads: 



^-E 
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Here, i, a and a are, respectively, the site (Mn or O), 
orbital (Mn- x^-y^, 3z^-l, O-p^;, Pj^, p^) and spin indices. 
The parameter £„ is the onsite energy of the orbital, 
(ij) indicates nearest neighbors, c^s are the creation op- 
erators, and Hi is the total number of electrons at i-th 
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FIG. 12: Orbitals considered in forming the Hamiltonian H|| 
appropriate for the out of the plane ferromagnetic configu- 
ration at the LMO/SMO interface. The parameter A is the 
onsite energy difference between the x^-y^ and 3z^-l orbitals. 
Depending on the strain condition A can be negative (com- 
pressive strain) or positive (tensile strain). The black dots 
indicate the occupied orbitals for one many-particle configu- 
ration, while the arrows indicate the spin states of the orbitals. 
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TABLE I: Slater-Koster tight binding hopping integrals be- 
tween the Mn-eg and O-p orbitals. In the three site model 
the value of Vpdcr is taken as 0.9 eV. 



site. The matrix elements Uajp are the Slater-Koster 
tight binding hopping integrals between the Mn-eg and 
0-p orbitals. We shall, for simplicity, take the Hund's 
coupling 3h as cxd so that only the eg states parallel to 
the t2g spin at a Mn site can be occupied. The symbol 
Si represents the t2g core spins and Sia is the spin of 
the Cg electron. The parameter 3sx represents the su- 
perexchange between the t2g core spins. Throughout this 
paper, we have taken Up — and Ud — 5 eV unless oth- 
erwise stated and also the on-site energy of the oxygen 
orbitals are taken as zero: ep — for all spins and all 
three p orbitals. 

The net exchange interaction between two Mn atoms is 
a sum of the antiferromagnetic superexchange Jsx, rnod- 
eled by the last term in the model Hamiltonian Eq. ([1]), 
and the double exchange 3dx mediated by the itinerant 
eg electrons, modeled by the remaining terms in the same 
equation, so that we have 



J = Jsx + Jdx ■ 



(2) 



The exchange interaction J is obtained by calculating the 
difference between the ground state energies correspond- 



ing to the ferromagnetic (FM) and the antiferromagnetic 
(AFM) alignment of the two t2g core spins: 
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Note that a positive (negative) value of J indicates an 
AFM (FM) interaction. Zsx is a simple additive term 
and it is, for the manganites, of the order of 26 meV^^. 
From the present model we will calculate the Jdx which 
depends on the occupancy of the eg states. 

Before we move on to the solution of the model, there 
is another point that needs to be made. The magnetic 
interaction between the planes, which is indicated by J in 
Fig. [l]may differ when we consider the exchange interac- 
tion between the Mn spins via the 0-p orbitals across the 
LaO plane or the SrO plane. However, from the charge 
contours {e.g. Fig. [6]), we see that there is very little dif- 
ference between the two oxygen atoms located on these 
planes, so that the Mn-O-Mn coupling may be expected 
to be nearly the same. So, in our model, we do not dif- 
ferentiate between these two interactions, so that J is the 
same across the LaO plane or the SrO plane. To distin- 
guish the exchange interaction out of the plane and in the 
plane, we have used the notation ^dx for the former and 
^'dx for the latter in the remaining part of this section. 



A. Out-of-plane exchange J 

First consider the out-of-plane exchange by evaluating 
the ground state energies for the ferromagnetic and the 



J 



antiferromagnetic configurations of the Mn-t2g spins of 
two Mn atoms along the z-axis. The model for the out- 
of-plane FM configuration is schematically shown in Fig. 
[T2] Listed in Table I are the Slater-Koster tight binding 
hopping integrals between Mn-e^ and 0-p orbitals, which 
shows that out of the three 0-p orbitals, only pz takes 
part in the hopping process along the z-axis. Therefore, 
in this model we have five spin- up orbitals (O-p^, two 
Mn-z-^-1 and two Mn- x^-y^) available for two spin-up 
electrons and one spin-down orbital (O-p^) available for 
the lone spin-down electron. The spin-down electron can 
only be on the O-pz orbital in the ferromagnetic case due 
to the infinite J/f and does not take part in the hopping 
process. Hence, we have a ten dimensional two particle 
configuration space i^C^ x^ Ci). 



We choose the two-particle (both particles with spins 
up) basis set in the order: \pdi), \d1d2), \d1d4), \pd2), 
\pdi), \d2di), \pd-i), \d2d3), \d3d4), and \d1d3), where p, 
di, d2, da and d4 respectively denote the 0-pz, Mnl-x^- 
y2, Mnl-3z2-l, Mn2-x2-y2, and Mn2-3z2-l orbitals. With 
this basis set, for the case of the two Mn-t2g core spins 
ferromagnetically aligned, the Hamiltonian for the itin- 
erant electrons H|| becomes: 
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Again, here Cd represents the onsite energy for the Mn- 
x^-y^ orbitals and A is the energy shift of the Mn-z-^-1 
orbital from the x^-y^ orbital due to strain. The onsite 
energy of the 0-p orbitals is taken as zero. Parameter t 
is the hopping matrix element {Ypda)- From DFT calcu- 
lations, we found that for compressive strain condition A 
is -l-ve and for tensile strain condition A is -ve. For the 
lattice matched interface the onsite energy of the x^-y^ 
and the 3z^-l orbitals are about the same. 

We now consider the antiferromagnetic case, where the 



two Mn-t2g spins are aligned antiferromagnetically. In 
this case, as seen from Fig. [T31 we see that there are six 
active orbitals, three spin-up and three spin-down. One 
can populate these orbitals with two spin-up electrons 
and one spin-down electrons or vice versa. (We do not 
consider the configurations where all three electrons have 
the same spins as this would correspond to an oxygen- 
to-manganese charge transfer state, which has a much 
higher energy.) In either of these cases we have a nine 
dimensional three particle configuration space {^C2 x'^ 
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Ci). These two sets of configurations do not interact 
with each other, as the model Hamiltonian Eq. ([1]) does 
not allow hopping between two opposite spins, and have 
the same ground-state energy. 
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FIG. 13: Orbitals considered in forming the Hamiltonian H|| 
appropriate for the out-of-plane exchange. Of the three oxy- 
gen p orbitals, only p^ has a non-zero hopping to the Mn eg 
orbitals. 



J 



By considering the configurations with two spin-up 
electrons and one spin-down electron and arranging the 
basis set in the order \ppdi), \pd1d2), \pd1d4), \d1d2d4), 
\pd1d3), \d1d2d3), \ppd2), \pd2d4), and \pd2d3), where the 
bar stands for the spin-down orbitals and unbar stands 
for the spin-up orbitals, the antiferromagnetic Hamilto- 
nian, H||, becomes: 



Hn - 



ea t t 

t 2ed + A + Ud t 

t 2ed + A t 

t t 3ed + 2A + Ud 
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ed + A t 

t 2ed + 2A 



.(5) 



[2ed + A] J 



r 



The above matrices can be diagonalized easily for all 
parameter values; however, for some limiting cases, one 
can solve these either analytically or using the pertur- 
bation theory, which then gives us considerable insights 
into the resulting exchange interactions. 

To this end, we first consider the limit when A is suf- 
ficiently large and negative. One immediately sees by in- 
specting the matrix that the ground state of iJ|| comes 
from the second block diagonal of Eq. (j4|) . Diagonaliza- 
tion then yields the ground state energy to be 



As seen from Eq. ([8]), the exchange energy is negative 
irrespective of the Hamiltonian parameters, indicating 
that the out-of-plane double exchange is ferromagnetic 
and robust. We have to add the superexchange term 
Jsx to this to get the net magnetic interaction. The re- 
sults are consistent with the numerically computed value 
of J presented in Fig. [TH which is discussed later. 

Similarly, we can also obtain an expression for a large 
and positive A. Again, it is immediately clear from the 
inspection of the Hamiltonians Eqs. (j4|) and ^ that the 

FM and AFM ground states come from the first sub-block 

Eij^ — 3(erf + A)/2 — V(e7+Ap + 8i^/2. (6) of H|| and H||. An analytical diagonahzation is not pos- 
sible in this case, however one can apply the fourth order 
For the AFM Hamiltonian i7|| in the same limit for A, non-degenerate perturbation theory to compute E^^ and 



the ground state comes from the third block diagonal of 
Eq. ([5]), yielding immediately the ground-state energy 



-Etj, the difference of which yields the result 
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The out-of-plane exchange energy 3ux is obtained by 
taking the difference between the FM and AFM ground 
state energies, so that for the case of the large and neg- 
ative values of A that we are considering, we get 

Jdx = Eij^ — E^i 



1/2 X ( V(q + A)2 + 4t2 - ^{ed + A)2 + 8*2). 

(8) 



We see that if the x2-y2 orbital is more occupied (i. e. A 
> 0), the out-of-plane magnetic ordering is AFM for any 
values of Cd and A. 

In Fig. [141 we have calculated the exchange interaction 
by simply a numerical diagonahzation of the two Hamil- 
tonian matrices for a general value of A, viz., -1.0 eV < 
A < 1.0 eV. The numerical results for J are consistent 
with with our analytical and perturbation results for the 
large values of A. When A is less than -0.15 eV we have 
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FIG. 14: Exchange parameter J {— Jdx + 3sx), with Jdx 
obtained from the numerical diagonalization of the Hamilto- 
nians (Eq. Q, ([Sjl, (|10|) . and dill) '), as a function of energy 
difference A between the x -y and the 3z -1 orbitals. For 
the calculation of J, the parameters td, U^, Vpdcr and Jsx are 
taken as 5 eV, 5 eV, 0.9 eV, and 26 meV respectively. When 
A is sufficiently positive (tensile strain) the magnetic inter- 
action is FM in the plane (J') and AFM between the planes 
(J). When A is sufficiently negative (compressive strain) the 
magnetic interaction reverses. 



a FM interaction and if it is greater than 0.15 eV we 
have an AFM interaction. For the intermediate region of 



A (-0.15 eV < A < 0.15 eV), the ground state swiftly 
changes between FM and AFM. 

Our model is consistent with the DFT calculations dis- 
cussed in the previous section (Fig. [3]), where we see that 
for compressive strain the 3z^-l orbitals are more occu- 
pied (A < 0) and the total energy calculation yields a 
out-of-plane FM configuration. If the strain is tensile, 
the x^-y^ orbital is more occupied (A > 0)and the out- 
of-plane magnetic ordering is AFM. 



B. In-plane exchange J' 

Now we consider the Mn-O-Mn exchange in the plane 
(xy-plane), i.e., between two Mn atoms on a Mn02 plane 
adjacent to the interface. In this case p^,, pj, and Cg are 
the only active orbitals since the hopping between p^ 
and eg orbitals is not allowed (Table I) . For concreteness 
we have taken the Mn-O-Mn to be along the x-axis and 
results are identical if it is along the y-axis. 

First consider the FM configuration of the two Mn- 
t2g spins. Here the itinerant electrons have a 10 dimen- 
sional configuration space as in the case of out-of-plane, 
but with Pj; as the active orbital instead of Pz. Taking 
the two-particle basis set in the order \pdi), \pd2), Iprfa), 
\pd4), \d1d2), \d1d3), \d1d4), \d2d3), \d2d4), and \d3d4) and 
denoting by t' and t" the hopping matrix elements 
and 



VsVpd. 



respectively, the Hamiltonian H|| becomes 
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Similar to out-of-plane AFM Hamiltonian, the in- 
plane AFM Hamiltonian has a nine dimensional con- 
figuration space. If we choose the basis set in the or- 



der \ppdi), \ppd2), \pd1d2), \pd1d3), \pd2d3), \d1d2d3), 
\pd1d4), \pd2d4), and \d1d2d4), the antiferromagnetic 
Hamiltonian H^i for the present case reads 
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TABLE 11: Mn-O-Mn exchange energy (in meV) in the plane 
(J') and out of the plane (J) obtained from the DFT calcula- 
tions and the three-site model for different values of the c/a 
ratio. The model uses the A values corresponding roughly to 
the DFT results for the three different strain conditions and 
the other parameters are: ^d = Ud ~ 5 eV, Vpda ~ 0.9 eV, 
and Jsx = 26 meV.-i^ 



An analytical diagonalization for both the ferromag- 
netic and antiferromagnetic Hamiltonian is not possible 
due to their non block diagonal nature and large dimen- 
sion. Non-degenerate perturbation theory can not be ap- 
plied as H|| contains degenerate states. However, in the 
limit A — > oo the ground state energies for the FM and 
AFM configuration is obtained from a 3 x 3 sub-matrix of 
H|l and a 2x2 sub-matrix of H|| respectively. The in- 
plane exchange energy J^x i^ the limit A ^ oo is then 
immediately obtained as 



JDX = W4 + '^t- 



iii2 



;2 + 8t"2)/2. (12) 



Here, we see that for large and positive value of A, inx 
is negative quantity which implies a FM interaction in 
the plane. 

Similarly, in the limit A ^ — oo, we find that 



J'dx = [^d + Ud + ^{€d + Udf+U"^]/2. 



(13) 



Eq. (|13p shows that for large and negative value of A, 
we have an AFM configuration in the plane. 

In Fig. [14] we have plotted the in-plane exchange J' 
as a function of A obtained from direct numerical diag- 
onalization of the full Hamiltonians Htt and Hti , which 



TT 
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indeed shows that for the in-plane exchange J', the mag- 
netic interaction in the plane switches from ferromagnetic 
to antiferromagnetic, as A is changed from positive to 
negative values. 

To summarize our analysis in the Section, we find that 
the strain-induced splitting in the Cg states is instrumen- 



tal in determining the magnetic properties at the inter- 
face. If x^-y^ orbital is relatively more occupied than the 
3z^-l orbital, the magnetic ordering at the LMO/SMO 
interface is more likely to be A-type with the FM and 
AFM configurations stabilized in the pane and out of the 
plane respectively. If the 3z^-l orbital is more occupied 
the magnetic ordering is more likely to be C-type with 
in-plane AFM ordering and out-of-plane FM ordering. 
When both the Cg orbitals are more or less equally occu- 
pied, the double exchange wins over the superexchange 
to stabilize the interface in a three dimensional FM con- 
figuration. 

The exchange interactions calculated from the three 
site model are in good agreement with the results ob- 
tained from DFT calculations as shown inTable II. 



V. EFFECT OF STRAIN ON A THICKER 
SUPERL ATTICE: {LaMnOa ) i / (SrMnOi ) 3 

In the preceding sections we have studied the ef- 
fects of strain on the magnetic interactions in the 
(LMO)i/(SMO)i superlattice. To generalize the strain 
effects on magnetism at the interface, in this section we 
have analyzed the magnetic properties of a thicker super- 
lattice, viz., (LM0)i/(SM0)3, which has both inner and 
interfacial Mn02 planes. 

In the (LM0)i/(SM0)3 superlattice, we have taken 
the G-type AFM configuration for the inner Mn02 planes 
as they belong to the SMO constituent of the superlat- 
tice. For the interfacial Mn02 planes, we have considered 
again the F-type, A-type and C-type magnetic configura- 
tions (Fig. [1]). In Fig. [121 we have shown the energetics 
for these three magnetic configurations as a function of 
the in-plane lattice parameter 'a'. From the figure we 
see that for lower values of the lattice constant 'a', the 
interface shows a stable C-type magnetic ordering and 
as we increase the value of 'a', the interface gradually 
moves towards an F-type magnetic configuration. For 
very high values of 'a', the interface stabilizes with the 
A-type magnetic configuration. 

For a strain free LMO/SMO superlattice, the in- 
plane lattice parameter coincides with the average lat- 
tice parameter acn (see section-II), which is 3.835 A for 
the (LM0)i/(SM0)3 superlattice. As in the case of 
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FIG. 15: Total energies for magnetic configurations A and 
C relative to the energy for the magnetic configuration F, 
as a function of the in-plane lattice parameter 'a' for the 
(LMO)i/(SMO)3 superlattice. The magnetic configurations 
A, F, and C are shown in Fig. [l] For strong tensile strain 
condition, we see a stable A-type magnetic ordering for the 
interfacial Mn02 layers, while for strong compressive strain 
condition the magnetic ordering is C-type. For the lattice 
matched case, the ferromagnetic ordering is stabilized in all 
directions. 



(LMO)i/(SMO)i superlattice (Fig. E]), here also we see 
that for strong compressive strain condition, the stable 
magnetic ordering at the interface is C-type, while for 
strong tensile strain condition the magnetic ordering is 
F-type. For the lattice matched interface we see a sta- 
ble FM ordering (F-type) in all directions. Our band 
structure calculation shows insulating behavior for strong 
compressive strain condition, while for lattice matched 
and tensile strain conditions it is metallic. 

To study the role of Mn-Cg orbitals on the interfacial 
magnetism in (LM0)i/(SM0)3 superlattice, in Fig. ^E\ 
we have shown the orbital projected Mn-Cg DOS for the 
interface Mn atoms for a = 3.92 A (tensile strain), 3.835 
A (lattice matched) and 3.75 A (compressive strain). 
From the figure we see that, analogous to the case of 
(LMO)i/(SMO)i superlattice, for the tensile strain con- 
dition the x^-y^ is relatively more occupied and 3z^-l 
orbital is less occupied, and opposite for the compres- 
sive strain condition. For the lattice matched interface, 
both the x^-y^ and 3z^-l orbitals are more or less equally 
occupied as before. 

Hence, as discussed earlier, for the tensile strain con- 
dition, we have a strong ferromagnetic double exchange 
coupling in the Mn02 plane through x^-y^ orbitals while 
between the planes we have a antiferromagnetic coupling 
due to superexchange between the t2g core spins. For 
the compressive strain condition the higher occupancy of 
Mn-3z^-l orbitals leads to a ferromagnetic coupling be- 
tween the planes and due to negligible occupancy of the 
x^-y^ orbital, the antiferromagnetic ordering is sustained 
in the Mn02 plane. For the lattice matched structure. 
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FIG. 16: Spin majority Mn-Cg DOS corresponding to the in- 
terface Mn atoms for the (LMO)i/(SMO)3 superlattice. The 
results are shown for the in-plane lattice parameter 'a' — 3.92 
A (tensile), 3.835 A (lattice matched) and 3.765 A (compres- 
sive). If the strain condition is tensile, the x^-y^ orbital is 
more occupied and the 3z^-l orbital is less occupied, and op- 
posite if the strain condition is compressive. For the lattice 
matched condition, both x'^-y'^ and 3z'^-l orbitals are more or 
less equally occupied. 



the orbital ordering is a combination of both x^-y^ and 
3z^-l orbitals and hence strong double exchange coupling 
both in the plane and out of the plane stabilizes the fer- 
romagnetic interaction in all directions. The similarity 
in the interfacial magnetic properties for both the super- 
lattices, viz., (LMO)i/(SMO)i and (LM0)i/(SM0)3 for 
different strain conditions suggests that the strain effect 
on magnetism at the LMO/SMO interface may be true 
for the general (LMO)„/(SMO)m superlattice as well. 



VI. SUMMARY 

In summary, we have studied the effect of the epi- 
taxial strain on the magnetic ordering at the interface 
of LMO/SMO superlattices by a detailed analysis on 
the (LMO)i/(SMO)i superlattice. We found that the 
epitaxial strain induces different orbital ordering which 
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in turn changes the magnetic ordering at the interface. 
The magnetic ordering at the interface is determined by 
the competition between the antiferromagnetic superex- 
change between the core t2g electrons and ferromagnetic 
double exchange between the itinerant e^ electrons. The 
strength of the latter in the Mn02 plane or between the 
planes strongly depends on the occupancy of the non- 
degenerate Bg orbitals. 

In case of a strong tensile strain condition, the higher 
occupancy of the x^-y^ orbital strengthens the double 
exchange coupling to stabilize the ferromagnetic order- 
ing in the Mn02 plane, while between the Mn02 planes 
the reduced double exchange coupling, due to the neg- 
ligibly occupied 3z^-l orbital, fails to overcome the an- 
tiferromagnetic t2g-t2g superexchange and stabilizes the 
A-type magnetic ordering at the interface. For strong 
compressive strain condition, the magnetic ordering re- 
verses, viz., higher occupancy of the 3z^-l orbital lead 
to a ferromagnetic coupling between the Mn02 planes. 



while the depleted x^-y^ occupancy allows the antifer- 
romagnetic ordering in the plane to make the C-type 
magnetic configuration as the most stable one. For a lat- 
tice matched structure, double exchange is strong enough 
both in the Mn02 plane and between the planes, due to 
more or less equally occupied x^-y^ and 3z^-l orbitals, to 
allow ferromagnetic ordering in all directions. 

The electronic structure calculations for the 
(LM0)i/(SM0)3 superlattice showed that the epi- 
taxial strain affects the magnetism at the interface in a 
similar way as the (LMO)i/(SMO)i superlattice. This 
suggests that the strain effect on magnetism may be 
similar for the general (LMO)„/(SMO)m superlattice. 
In addition, similar considerations regarding the effects 
of strain on orbital ordering and magnetism should be 
valid for interfaces between other perovskite oxides as 
well. 
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